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I . The equilibrium phase diagram of a two-dimensional Ising model with competing exchange and 

■ dipolar interactions is analyzed using a Monte Carlo simulation technique. We consider the low 
temperature region of the (5, T) phase diagram (5 being the ratio between the strengths of the 

£N| ■ exchange and dipolar interactions) for the range of values of S where striped phases with widths 

' h — 1 and h = 2 are present. We show that the transition line between both phases is a first order 

one. We also show that, associated with the hrst order phase transition, there appear metastable 
states of the phase ft = 2 in the region where the phase ft = 1 is the thermodynamically stable one 
' and viceversa. 

o\ ■ 

! ,! PACS numbers: 75.40. Gb, 75.40.Mg, 75.10.Hk 

The physics of ultrathin films and quasi-two-dimensional systems is of current interest because of its many tech- 
nological applications. For instance, metal-on-metal films are used in electronics, data storage and catalysis. In 
. » ■ particular, the use of ultrathin magnetic films as data storage devices requires a high degree of accuracy and spatial 
resolution in the magnetization control. Hence, a lot of experimental and theoretical effort has been devoted in the last 

■ yearsu to understand both the equlibrium and out-of-equilibrium properties of these kind of systems. Realistic theo- 
Jh [ retical descriptions of these systems include both the exchange and the dipolar interactions between the microscopic 

spins. The competition between long ranged antiferromagnetic dipolar interactions and short ranged ferromagnetic 
exchange interactions can give rise to a variety of unusual and interesting macroscopic phenomena. Works in two 
dimensional uniaxial spin systems, where the spins are oriented perpendicular to the lattice and coupled with these 
kind of interactions, have shown a very rich phenomenological scenario concerning both its .equilibrium statistical 
O \ mechanicsou and non-equilibrium dynamical propertiesoU. In particular, some of these resultsEl showed the existence 
of different types of slow relaxation dynamics when the system is quenched from a high-temperature configuration to 
a subcritical temperature, depending on the relative strenghts of the dipolar and exchange interactions. This change 
in the relaxation properties could be an effect of a change in the low temperature equilibrium properties of the system, 
such as the presence of metastable states. In this work we investigate in detail the low temperature phase diagram of 
, that system in the region where the change in the relaxation properties has been observed. 
The above mentioned system is described by the Ising like Hamiltonian 

where the spin variable <Ji = ±1 is located at the site i of a square lattice, the sum X]<ij> runs over all pairs of 
nearest neighbor sites and the sum y\. ~ runs over all distinct pair of sites of the lattice; is the distance (in crystal 
units) between sites i and j; 6 represents the quotient between the exchange Jo and dipolar Jd coupling parameters, 
where the energy is measured in units of Jd, which is always assumed to be antiferromagnetic (Jd > 0). Hence, 6 > 
q | means ferromagnetic exchange coupling. 

The overall features of the the finite temperature phase diagram of this model were described by Maclsaac and 
coauthorsa by means of Monte Carlo calculations on 16 x 16 lattices and analytic calculations of the ground stateQ. 
They found that the ground state of Hamiltonian ([!]) is the antiferromagnetic state for 5 < 0.425. For S > 0.425 
the antiferromagnetic state becomes unstable with respect to the formation of striped domain structures, that is, 
to state configurations with spins aligned along a particular axis forming a ferromagnetic stripe of constant width 
ft, so that spins in adjacent stripes are anti-aligned, forming a super lattice in the direction perpendicular to the 
stripes. They also showed that striped states of increasingly higher thickness ft become more stable as 5 increases 
from 5 = 0.425. Moreover, they showed that the striped states are also more stable than the ferromagnetic one for 
arbitrary large values of S, suggesting such a phaseito be the ground state of the model for S > 0.425. Monte Carlo 
calculations on finite lattices at low temperature&Q gave further support to this proposal, at least for intermediate 
values of S. Furthermore, such simulations have shown that striped phases of increasingly higher values of ft may 
become thermodynamically stable at finite temperatures for intermediate values of 6. These results are in agreement 
with other analytic onescm. For low values of S the system presents an antiferromagnetic phase at low temperatures. 
At high temperatures, of course, the system always becomes paramagnetic. Specific heat calculations showed that the 
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transition between the paramagnetic and the striped phases is a second order oneH, while the nature of the transitions 
between the different striped phases was not clearly determined. 

We performed Monte Carlo simulations of Hamiltonian (]l|) on square lattices up to 48 x 48 sites using periodic 
boundary conditions and heat bath dynamics. Our calculations focused on the low temperature region of the (S,T) 
phase diagram for values of 5 between 0.4 and 2, which includes the transition line between the h — 1 (hi) and the 
h = 2 (h2) striped phases. 

First we calculated through the energy fluctuations the specific heat C as a function of the temperature for different 
values of S and different system sizes up to 48 x 48 sites. The typical behaviour of C is shown in Fig|l] for 5 = 1.1. 
By considering the peaks in the specific heat we obtained the second order critical line between the paramagnetic and 
the low temperature ordered phases (hi and h2). These results (triangles in Fig.[|) slightly improve those obtained by 
Maclsaac and coauthorsa for 16 x 16 lattices, thus showing a fast convergence of the critical temperature for increasing 
sistem sizes, at least for small values of 6. 
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FIG. 1. Specific heat C vs. temperature T for 8 — 1.1 and different system sizes N. 



Next, we analized the transition between the hi and h2 phases. To do this we introduced the staggered magnetization 

1 N 

m hi = jj ( 2 ) 

x,y=l 

where N is the total number of lattice sites; rrihi equals one in the ground state corresponding to the striped phase 
with h — 1. We calculated the thermal average Mhi =< m/ii > as a function of T as well as the fluctuations of this 
quantity through the associated staggered susceptibility 

Xhi = {m 2 hl ) - (m hl ) 2 . (3) 




T T 



2 



FIG. 2. Staggered magnetization Mhi {Mh2) an d associated susceptibility Xhi (xu) vs. T when the system is heated starting 
at T = from an initial configuration in the corresponding ground state hi (h2), for N — 24 x 24. (a) 8 = 1 corresponding to 
the stable phase hi; (b) 8 = 2 corresponding to the stable phase h2. 

In a similar way we defined the staggered magnetization Mhi and its associated susceptibility Xh2, corresponding to 
the striped phase h2. We simulated the heating of the system from T = to a temperature higher than the critical 
one for different values of 8 and starting from different initial configurations. To analize the stability of both phases 
in the different parts of the phase diagram we first calculated the staggered magnetization Mhi and susceptibility Xhi 
starting from a striped configuration hi; we also calculated the staggered magnetization Mh2 and susceptibility Xh2 
starting in a striped configuration h2. The typical behaviours are shown in Fig. 2 for 6=1 (Fig. 2a), where the ground 
state corresponds to hi, and for 8 = 2 (Fig. 2b), where the ground state corresponds to h2. For 8 = 1 we see that Mhi 
falls down to zero at the critical temperature while Xhi shows a sharp peak at the same temperature, as expected 
for the order parameter and its conjugated susceptibility at this type of transition. On the other hand, Mh2 remains 
stable up to some temperature T2 < T c , where it losses stability and falls down abrubtly to zero, while Xh2 is zero 
for almost every temperature, except near T = T2. A similar behavior is observed for 8 = 2, but with (Mhi,Xhi) and 
{Mh2,Xh2) interchanged. In this case the phase hi losses stability at some temperature Ti < T C (S). 

All these results show the existence at low temperatures of metastable states h2 in the region of the phase diagram 
corresponding to the stable phase hi and viceversa. The numerical calculation of stability lines Ti(8) and T2(S) is 
depicted in Fig.[|by means of circles and squares respectively. We see that these lines cross at a certain value S ~ 1.26 
where they join smoothly with the second order critical lines. 

All these results suggest the existence of a first order phase transition line separating the phases hi and h2. To 
verify this we considered the free energy 
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FIG. 3. Free energy F of the phases hi and h2 vs. 8 for T = 0.2. 

where U = (H) . We then calculated the free energy of both phases hi and h2 by heating the system from T = up 
to a fixed value of temperature T, starting from two different initial configurations corresponding to the ground state 
of hi and h2, and for different values of 8. In Fig. |] we see an example of this calculation for T = 0.2. The free energy 
of the phase hi was calculated for increasing values of S, starting from some small value well inside the region where 
this phase is stable, up to the the value of 8 corresponding to the stability line of hi at the given temperature. The 
same calculation was repeated a the same temperature for the free energy of the phase h2, but for values of 8 ranging 
from the stability line up to some value well inside the region where h2 is stable. The observed continuous change of 
the minimal free energy from one phase to the other, with a discontinuous change in the slope as 8 is varied is a clear 
evidence of a first order phase transition. Also the multivalued nature of the free energy gives further evidence of 
the metastable nature of these phases in some parts of the phase diagram. Repeating these calculations for different 
values of T we obtained the almost vertical first order transition line between the two phases, showed in Fig.|] by 
means of diamonds. The shaded region in Fig.H indicates the presence of metastable states. 
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FIG. 4. Numerical calculation of the phase diagram (5, T) in the region of parameters under study. Filled triangles correspond 
to the critical temperatures T c (5) obtained by specific heat calculations for the phase transition between the ordered striped 
phases hi and h2 and the paramagnetic one. Filled circles (open squares) correspond to the stability line of the hi (h2) phase, 
obtained by analizing the staggered magnetization Mux {Mki}- Filled diamonds correspond to the first order transition line 
between the hi and h2 phases, obtained by free energy numerical calculations. The shaded region correspond to the presence 
of metastable states. 



Finally, we considered the relaxation of the system in the metastable region starting from a non-uniform initial 
configuration. We prepared the system in a configuration with one half of the system in the hi phase and the other 
half in the h2 phase, as shown in FigM This particular configuration facilitates the nucleation of the stable phase 
whatever it be (hi or h2) . In the same figure we show also the time evolution of the energy for two different values of 
5 located at both sides of the transition line S = 1.26. We see that the energy evolves towards the mean energy of hi 
{Ehi) at the corresponding temperature when S < 1.26, while it relaxes to the mean energy of h2 (£^2) when S > 1.26. 
This shows that indeed the metastable configurations are unstable agains the nucleation of the corresponding stable 
phase. 
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FIG. 5. (a) Initial configuration with 50% in the ground state corresponding to hi and 50% in the ground state corresponding 
to h2. (b) Time evolution of the energy of the system E starting from the initial configuration of figure (a), for T = 0.2 and 
N = 24 x 24. 



Summarizing, we presented numerical evidence of the existence of metastable states in the model described by 
Hamiltonian (|lj) in the low temperature region of the (6, T) phase diagram where the transition between hi and h2 
phases takes place. We also showed that these metastable states are associated with a first order phase transition line 
between both phases. The presence of these metastable states may alter the normal domain growth process when the 
system is quenched from high temperatures into the ordered phase, by pinning the walls of the stable phase domains. 
Such process could slow down the normal coarsening dynamics, depending on whether the quench drops the system 
into the metastable region or not, thus explaining the observed change in the relaxation dynamics in this region of the 
phase spaced. Some work along this line is in progress and will be published elsewhere. Another point that deserves 
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further investigation is the possible presence of more complex metastable states for higher values of 5, as observed 
in a related three-dimensional model with [Competing nearest neighbors ferromagnetic interactions and long-range 
antiferromagnetic Coulomb-like interactionsu. 
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